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Abstract 

The continuum percolation for Markov (or Gibbs) germ-grain models is investigated. 
The grains are assumed circular with random radii on a compact support. The morpho- 
logical interaction is the so-called Quermass interaction defined by a linear combination 
of the classical Minkowski functionals (area, perimeter and Euler-Poincare characteristic). 
We show that the percolation occurs for any coefficient of this linear combination and for a 
large enough activity parameter. An application to the phase transition of the multi-type 
Quermass model is given. 

KEY-WORDS: Stochastic geometry, Gibbs point process, germ-grain model, Quermass 
interaction, percolation, phase transition. 
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1 Introduction 



The germ-grain model is built by unifying random convex sets- the grains -centered at 
the points- the germs -of a spatial point process. It is used for modelling random surfaces 
and interfaces, geometrical structures growing from germs, etc. For such models, the con- 
tinuum percolation refers mainly to the existence of an unbounded connected component. 
This phenomenon expresses some macroscopic properties of materials as permeability, con- 
ductivity, etc. Moreover, it turns out to be an efficient tool to exhibit phase transition 
in Statistical Mechanics [21 0]. For theses reasons, the continuum percolation has been 
abundantly studied since the eighties and the pioneer paper of Hall [8|. 

When the grains are independent and identically distributed, and the germs are given 
by the locations of a Poisson point process (PPP), the germ-grain model is known as the 
Boolean model. In this context, the continuum percolation is well understood; see the 
book of Meester and Roy |13| for a very complete reference. One of the first results is the 
existence of a percolation threshold z* for the intensity parameter z of the stationary PPP: 
provided the mean volume of the grain is finite, percolation occurs for z > z* and not for 
z < z . 

Because of the independence properties of the PPP, the Boolean model is sometimes 
caricatural for the applications in Biology or Physics. Mecke and its coauthors [111112] have 
mentioned the need of developing, via Markov or Gibbs process, an interacting germ-grain 
model in which the interaction would locally depend on the geometry of the set. For this 
purpose, let us cite the Widom-Rowlinson model |16) . the area interaction process [T| and 
the morphological model |12j . Thus, Kendall, Van Lieshout and Baddeley suggested in 
[9] a generalization of the previous models, called the Quermass Interaction Process. In 
this model, the formal Hamiltonian is a linear combination of the fundamental Minkowski 
functionals, namely in M. 2 the area A, the perimeter C and the Euler-Poincare characteristic 
X- 

h = o 1 a + e 2 c + e 3X . 

The existence of infinite volume Gibbs point processes for the Hamiltonian H has been 
recently proved in [3]. This paper focuses on the continuum percolation for such processes. 

The existence of a percolation threshold z* for the Boolean model relies on a basic 
(but essential) monotonicity argument: see [13] , Chapter 2.2. This argument fails in the 
case of Gibbs point processes with Hamiltonian H. So, no percolation threshold can be 
expected in our context. However, other stochastic arguments as stochastic domination or 
FKG lead to percolation results. In [2], Chayes et al prove that percolation occurs for z 
large enough and 62 = 63 = 0. To our knowledge, the percolation phenomenon for other 
values of parameters 61,62,63 has not been investigated yet. 

Our main result (Theorem [1]) states that, for any 61,62, 63 (positive or negative), per- 
colation occurs with probability 1 for z large enough. The only assumption bears on the 
random radii of the circular grains: they have to belong to a compact set not containing 
0. The proof of this theorem is relatively easy in the case 63 = 0. Indeed, the local energy 
h((x, R),uj)- the energy variation when the grain B{x, R) is added to the configuration oj - 
is uniformly bounded (see Lemma [4.12p and by classical stochastic comparison with respect 
to the Poisson Process the result follows. So the main challenge of the present paper is the 
proof of Theorem Q] when 63 7^ 0. In this setting, the local energy becomes unbounded from 
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above and below and classical stochastic comparison arguments for point processes fail. In 
interpreting the percolation in our model via a site percolation model (see the beginning 
of Section [4]) , we prove the result thanks to a stochastic domination result for graphs due 
to Liggett et al (Theorem 1.3 in |10| or Lemma 14.21 below). An arduous control of the 
hole number variation, when a new grain is added, is the main technical issue. We prove 
essentially that this variation is moderate for a large enough set of admissible locations of 
grains. Let us mention that our proof is inspired by the one of Proposition 3.1 in pQ. 

Following [2j[3], we use our percolation result (Theorem [1]) to exhibit a phase transition 
phenomenon for Quermass model with several type of particles (Theorem [2]) . 

Our paper is organized as follows. In Section [2j the Quermass model and the main 
notations are introduced. The local energy h((x, R),u>) is defined in (|2.3|) . Section [3] 
contains the results of the paper. Section T3. 21 is devoted to the case 63 = and Section I3TB1 
to the phase transition result. The proof of Theorem [1] is developed in Section [H 

2 Quermass Model 

2.1 Notations 

We denote by B(M. 2 ) the set of bounded Borel sets in R 2 with a positive Lebesgue measure. 
For any A and A in B(M?), A © A stands for the Minkoswki sum of these sets. Let 
< Ro < R\ be some positive reals and £ be the product space IR 2 x [Rq,Ri] endowed 
with its natural Euclidean Borel cr-algebra cr(£). For any A £ £>(IR 2 ), £\ denotes the space 
A x [Rq,Ri]. A configuration a; is a subset of £ which is locally finite with respect to its 
first coordinate: ^(ujn£\) is finite for any A in B(M?). The configuration set f] is endowed 
with the c-algebra T generated by the functions oj 1— >■ #(w n A) for any A in a(£ ). 
We will merely denote by uj\ instead of u Pi £\ the restriction of the configuration uj (with 
respect to its first coordinate) to A. Moreover, for any (2, R) in £, we will write u) U (cc, R) 
instead of co U {(x, R)}- 

A configuration uj S f2 can be interpreted as a marked configuration on BL 2 with marks 
in [Ro,R\]. To each (x,R) € a; is associated the closed ball B(x,R) (the grain) centered 
at x (the germ) with radius R. The germ-grain surface u> is defined as 

lo= (J B(x,R) . 

2.2 Quermass interaction 

Let us define the Quermass interaction as in Kendall et al. |9] . The energy (or Hamiltonian) 
of a finite configuration u) in Q is defined by 

H(u) = 0i.A(w) + 9 2 C(Q) + 6 3X (u) , (2.1) 

where 9\, 62 and 83 are three real numbers, and A, C and x are the three fundamen- 
tal Minkowski functionals, respectively area, perimeter and Euler-Poincare characteristic. 
This last one is the difference between the number of connected components and the num- 
ber of holes. Recall that a hole of Q is a bounded connected component of u) c . Hadwiger's 
Theorem ensures that any functional F defined on the space of finite unions of convex 
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compact sets, which is continuous for the Hausdorff topology, invariant under isometric 
transformations and additive (i.e. F(A U B) = F(A) + F(B) - F(A n B)) can be de- 
composed as in (|2.ip . This universal representation justifies the choice of the Quermass 
interaction for modelling mesoscopic random surfaces fl"2] . 

The energy inside A £ B(M. 2 ) of any given configuration co in (finite or not) is defined 

by 

H a (oj) = H(cj a ) - H(cj A \ A ) , (2.2) 

where A is any subset of M? containing A © B(0, 2R\). By additivity of functionals A, C 
and \, the difference H a (oj) does not depend on the chosen set A. 

Let us end with defining the local energy h{(x,R),uj) of the marked point (x,R) £ E 
(or of the associated ball B{x,R)) with respect to the configuration co: 

h((x,R),uj) = H a (ujU(x,R)) - H a (uj) , (2.3) 

for any A £ B(M?) containing x. Remark this definition does not depend on the choice 
of the set A. The local energy h((x,R),u) represents the energy variation when the ball 
B(x, R) is added to the configuration u. 



2.3 The Gibbs property 

Let Q be a reference probability measure on [Rq,Ri]. Without loss of generality, Rq and 
Ri can be chosen such that, for every e > 0, 

Q([Ro + e, R{\) < 1 and Q{[R , R x - e\) < 1 . (2.4) 

Let z > 0. Let us denote by A the Lebesgue measure on M 2 and by 7r z the PPP on 
E with intensity measure zX <S> Q- Under 7r 2 , the law of the random surface ui is the 
stationary boolean model with intensity z > and distribution of radius Q. Finally, for 
any A £ B(M. 2 ), let us denote by tt a the PPP on £ A with intensity measure z\ A <giQ, where 
Aa is the restriction of the Lebesgue measure A to A. 

Definition 2.1. A probability measure P on £1 is a Quermass Process for the intensity 
z > and the parameters 61,62, #3 if P is stationary and if for every A in B(M?), for every 
bounded positive measurable function f from f2 to R ; 

f(u)P(du;) = J J / Ku WA c)^-^e-^Ku- A c)^ ((iLi; / v) p (da;) ^ (25) 

where Z A (u A c) is the partition function 

Z A (u Aa ) =\l e-^K^) 4 (^) . 



The equations (|2.5j) for all A £ B(M?) -are called DLR for Dobrushin, Landford and 
Ruelle. They are equivalent to: for any A £ B(M. 2 ), the law of uj a under P given uj A c is 
absolutely continuous with respect to the Poisson Process 7r A with the local density 

g A {J A \uj A c) = ^ye-^Hu^) . (2. 6) 
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See |15| for a general presentation of Gibbs measures and DLR equations. 

The existence, the uniqueness or non-uniqueness (phase transition) of Quermass pro- 
cesses are difficult problems in statistical mechanics. The existence has been proved re- 
cently in [3], Theorem 2.1 for any parameters z > and 61,82, 63 in R . A phase transition 
result is proved in [2 [6j [16] for Rq = R±, 62 = 63 = and for 6\ = z large enough. 

3 Results 

3.1 Percolation occurs 

We say that percolation occurs for a given configuration cj 6 Q if the subset u of M 2 con- 
tains at least one unbounded connected component. The set of configurations such that 
percolation occurs is a translation invariant event. Its probability, called the percolation 
probability, equals to or 1 for any ergodic Quermass process. However, the Quermass 
processes are not necessarily ergodic (they are only stationary) and their percolation proba- 
bilities may be different from and 1. Besides, in [2], Chayes et al have built two Quermass 
processes, both corresponding to 62 = 63 = and 8\= z large enough, whose percolation 
probabilities respectively equal to and 1. Since any mixture of these two processes is still 
a Quermass process, the authors obtain Quermass processes whose percolation probabili- 
ties equal to any value between and 1. 

Our main result states that percolation occurs with probability 1 for any (ergodic or not) 
Quermass process whenever the intensity z is large enough. 

Theorem 1. Let Rq > and 61,62,0-3 6 R. There exists z* > such that for any 
Quermass process P associated to the parameters 61,62, 63 and z > z* , percolation occurs 
P-almost surely. 

The proof of Theorem [T] is based on a discretization argument which allows to reduce 
the percolation problem from the (continuum) Quermass model to a site percolation model 
on the lattice 1? (up to a scale factor). This proof is rather long and technical so it is 
addressed in Section [H 

Let us point out here that our theorem does not claim z* is a percolation threshold. In 
other words, for z < z* , the percolation may be lost and recovered on different successive 
ranges. 

Another natural question involves the number of unbounded connected components. Fol- 
lowing the classical arguments for continuum percolation, we prove that this number is 
almost surely equal to zero or one. 

Proposition 3.1. For any Quermass process P the number of unbounded connected com- 
ponent is a random variable in {0, 1}. 

Proof. It is well-known that any Gibbs measure is a mixture of extremal ergodic Gibbs 
measures. For each ergodic Quermass process P, the number of connected component is 
almost surely a constant in N U {+00}. For any A € £?(M 2 ), thanks to the DLR equations 
(12. 5p . it is easy to prove that the law of u>\ under P is equivalent to 7r^. Therefore, in 
following the general scheme of the proof of Theorem 2.1 in |13| . we show that the number 
of connected components is necessary or 1. □ 
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3.2 Percolation when #3 = 

In the particular case 83 = 0, Theorem [1] can be completed and proved in a simple way. 

First, let us recall the definitions involving the stochastic domination for point pro- 
cesses. We follow the notations given in [5]. An event A in T is called increasing if for 
every uj G A and any uj' G 0, containing uj then uj' £ A too. Let P and P' be two proba- 
bility measures on Vt. We say that P is dominated by P' , denoted by P ^ P' , if for every 
increasing event A £ J 7 , P{A) < P'{A). In this section, we focus our attention on the 
increasing event "there exists an unbounded connected component". 

Let P be any Quermass process and assume #3 = and Rq > 0. Thanks to Lemma 
I4.12[ the local energy can be uniformly bound: there exist constants Co and C\ such that 
for any (x, R) G £ and uj G ft, 

C < h((x,R),w) < Ci . (3.1) 
Combining (|3.ip and Theorem 1.1 in [5], we get the following stochastic dominations: 

ze~ c i j T3 j ^ze~ c 

Now, the (stationary) Boolean models corresponding to -zr 26 Cl and 7r ze °° admit positive 
and finite percolation thresholds (see |14| . Chapter 3). It follows : 

Proposition 3.2. Let Rq > 0. For every 81,62 in K, there exist constants zq,z\ such that 
for any Quermas Process P associated to parameters z, 81,82 and 83 = 0, the percolation 
occurs P -almost surely if z> Z\ and does not occur P -almost surely if z < zq. 

Proposition 13.21 improves Theorem [1] in the case 83 = since it ensures the existence 
of a subcritical regime. 

It is worth pointing out here that the uniform bounds in (|3.ip do not hold whenever 
83 7^ 0. Precisely, this is the hole number variation which cannot be uniformly bounded. 

3.3 Phase transition for multi-type Quermass Process 

In this section, the multi-type Quermass model is introduced and a phase transition is 
exhibited, i.e. the existence of several Gibbs processes for the same parameters is proved. 

Let K be a positive integer. The X-type Quermass model is defined on the space Qk 
of configurations in £k = M 2 x [Rq, R\] x {1,2, . . . ,K}. Each disc is now marked by a 
number specifying its type. We don't give the natural extension of the notations involving 
the sigma-field and so on. 

The following Quermass energy function is defined such that all discs of a connected com- 
ponent have the same number. This is a non-overlapping multi-type germ-grain model. 
Precisely the energy of a finite configuration uj is now given by 

H(uj) = 0iA(u) + 8 2 £{uj) + 8 3X (uj) + 4>{\x-y\-R-R') , (3.2) 

(x,R,i),(y,R',j)£u> 

where <f> is an hardcore potential equals to infinity on ] — 00, 0] and zero on ]0, +00]. The 
energy inside A G B(M?) of any finite or infinite configuration uj is defined as in (|2.2p with 
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the convention +00 — 00 = +00. The definition of the i^-type Quermass process via the 
DLR equations follows as in Definition 12.11 

The proof of the existence of such processes is similar to the one of the existence of 
Quermass process. See Theorem 2.1 of [3] for more details. Here is our phase transition 
result: 

Theorem 2. Let Rq > 0. For any Q\, 62 and #3 in M, there exists zq > such that, for 
any z > zq, there exist several K-type Quermass Processes. The phase transition occurs. 

The proof essentially follows the scheme of the one of Theorem 2.2 of [2] or Theorem 1.1 
of [1]. It is based on a random-cluster representation (or Gray Representation) analogous to 
the Fortuin-Kasteleyn representation of the Potts model. The existence of an unbounded 
connected component allows to prove the existence of a K-type Quermass process in which 
the density of particles of a given type is larger than the ones of the other types. By 
symmetry of the types, we prove the existence of at least K different if- type Quermass 
processes. 

4 Proof of Theorem Q] 
4.1 General scheme 

In the following, P denotes a stationary Quermass process on Q associated to the intensity 
z > and the parameters 6\, 62, #3 G R. 

Let £ be a real number such that £ > 2R\ + 2Rq. Let us define the diamond box A as 
the interior of the convex hull of the eight points (3£,0), (6£,0), (9£,3£), (9£,6£), (6£,9£), 
(3£,9£), (0,6^) and (0,3£). This large octagon contains four smaller boxes B N , Bs, B E 
and By with side length £; precisely B N = (A£,7£) + [0,£] 2 , B s = (4£,£) + [0,£] 2 , B E = 
(7£, U) + [0, £f and B v = {£, U) + [0, £] 2 . The subscripts N, S, E and W refer to the cardinal 
directions. See Figure [TJ Thus, let us introduce the indicator function £ defined on and 
equal to 1 if and only if the two following conditions are satisfied: 

(CI) Each box i? N , Bs, B E and By, contains at least one point of wa; 

(C2) The number N^(oj) of connected components of (Da having at least one ball centered 
in one of the boxes B^, Bs, B E or B v , is equal to 1. 

In other words, £(w) = 1 means the boxes B®, Bs, B E and B v are connected through Da- 
For any x G (6£Z) 2 , let t x be the translation operator on the configuration set £ defined 
by {y, R) G t x uj if and only if (y + x,R) G 00. Hence, we can define the translated indicator 
function £ x of £ on the translated box A x = x + A by £ x -(cl>) = £(t x lu). Let us remark that 
£ x (w) only depends on the restriction of the configuration oj to the box A x . Moreover, 
thanks to the stationary character of the Quermass process P, the random variables ^ x , 
x G (Q£7 J ) 2 , are identically distributed. They are dependent too. 

Let us consider x,y G (6£Z) 2 such that y = (6£, 0) + x. The boxes A x and A^ have in 
common a cardinal box, i.e. x + B E = y + B v . So, the condition £, x (u>) = £^(w) = 1 ensures 
that the cardinal boxes of A^ and A y are connected together through the restriction of ui 
to A x U Ay. The same is true when y = (0, 6£) + x. This induces a graph structure on the 
vertex set V = (6£'Z) 2 : for any x, y G V , {x, y} belongs to the edge set E if and only if 

y-xe {±(6£,0),±(0,6£)} . 
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U 6£ 9£ 



Figure 1: Here is the diamond box A. The light gray set represents the configuration u restricted 
to A. The dark gray squares are the fourth cardinal boxes B N , B s , B E and B v with side length £. 
On this picture, conditions (CI) and (C2) are fulfilled, i.e. = 1. 

The graph (V, E) is merely the square lattice Z 2 with the scale factor 6£. The family 
{£x,% G V} provides a site percolation process on the graph (V,E). It has been built so 
as to satisfy the following statement. 

Lemma 4.1. Let u) £ £1 such that percolation occurs in the site percolation process {£ x , x S 
V}. Then, so does for uj. 

Let Hp be the Bernoulli (with parameter p) product measure on {0, 1}^. A stochastic 
domination result of Liggett et al [10J (Theorem 1.3) allows to compare the site percolation 
processes induced by the family {£ x ,% £ V} an d n p . Here is an adaptated version to our 
context. Basic definitions about stochastic domination for lattice state spaces are not recall 
here. They are similar to the ones presented in Section 13.21 for point processes. See also 

Lemma 4.2. Let p £ [0, 1]. Assume that, for any vertex x G V , 

P{t x = l\i y :{x,y} £E)>p a.s. (4.1) 

Then the distribution of the family {£ X y x £ ^} stochastically dominates the probability 
measure n^( p ), where f : [0, 1] — > [0, 1] is a deterministic function such that f(p) tends to 
1 as p tends to 1 . 

Actually, Theorem [1] straight derives from Lemmas 14.11 and 14.21 Let us first recall that 
in the site percolation model on the graph (V, E), there exists a threshold value p* < 1 such 
that percolation occurs with LLp-probability 1 whenever p > p* . See the book of Grimmett 
[7], p. 25. So, let p be a real number in [0, 1] such that 

f(p) > p* . (4.2) 



8 



Whenever the Quermass process P satisfies (|4.ip for that p, then combining Lemmas 14.11 
and 14.21 percolation occurs P-a.s. Therefore it remains to show that for any p > 0, hypoth- 
esis (|4.ip holds for z large enough. 

The next result claims that each Borel set of R 2 , sufficiently thick in some sense, 
contains at least one element of the configuration u with a probability tending to 1 as the 
intensity z tends to infinity. It will be proved at the end of this section. 

Lemma 4.3. Let l^cM 2 such that there exist U G B(M. 2 ) with positive Lebesgue measure 
and e > satisfying U © -6(0, R\ + Rq + e) C V . Then there exists a constant C > 0, 
depending on X(U) and e, such that for any configuration uj G £1 and for any z > ; 

P(u v =0|w V c) < Cz^ 1 . 

Since the Quermass process P is stationary, it is sufficient to prove (|4.ip with x = (0, 0). 
So, we focus our attention on the diamond box A = A( 0) o) an d use Lemma 14.31 to check 
that condition (CI) is fulfilled in this box. Since i? n , Bs, B E and By are sufficiently thick 
(with side length £ > 2R\ + 2R ), it follows 

P {uj Bi =0|w A c) = P(P (u Bi = | cJBf) I w Ac ) < Cz" 1 , 

for any % G {N,S,E, W}. So the conditional probability that oj satisfies (CI) is larger than 
1 -ACz~ l . 

The equation N^.(uj) = forces the box i? N (for instance) to be empty of points of the 
configuration uj. Hence, 

P{N*{ui) = 0|wa=) < Cz^ 1 . 

Checking that condition (C2) is fulfilled in the diamond box A needs what we call the Con- 
nection Lemma (Lemma 14. 4p . This result states the conditional probability that N^,(lj) 
is larger than 2 converges to uniformly on the configuration outside A. This is the heart 
of the proof of Theorem [TJ Its technical proof is given in Section 14.21 

Lemma 4.4 (The Connection Lemma). There exists a constant C > such that for any 
configuration uj G and for any z > 0, 

P (N^{lj) > 2 | u A c) < C'z- 1 . (4.3) 

The above inequalities and the Connection Lemma imply that conditions (CI) and 
(C2) are fulfilled in A with a probability tending to 1 as z tends to oo: 

P (foo) M = 1 1 oja°) > 1 - (5C + C')z- 1 . 

The hypothesis (|4.ip then follows. Let x be a vertex of the graph (V, E) which is not a 
neighbor of (0,0). By construction, the box A x is included in A c = A^ 0Q ^ (since A is an 
open set). This means the random variable £ x is measurable with respect to the u-algebra 
induced by the configurations restricted to A^ Q( ^. So, 

P (£(o,o) = 1 1 : {(0, 0), x} £E) > 1 - (5C + C')z~ x , 

and the hypothesis (|4.ip holds with x = (0,0) and any p G [0, 1[, provided the intensity z 
is large enough. This ends the proof of Theorem [TJ 
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Proof. (Lemma l4.3p Let U £ B(M. 2 ) be a bounded Borel set with positive Lebesgue measure 
and V D U B(0, R 1 + R + e). First, let us write: 

P(av = 0|o>H = 7 } J K v=$ e- H y^y^^(duj v ) 
Z, v [w V c) J Qy 

e -z\(V) 

(4.4) 



(wy c ) 

since the empty configuration has a null energy, i.e. Hy(ooyc) = 0. A configuration to 
whose restriction to V satisfies 

#^Ux[Ro,Ro+e] = 1 and ^V\U = 

is reduced to a ball B(x, R) centered at a x in U and with a radius Rq < R < Rq + e. Since 
the ball B{x,R) does not overlap u)ye, its energy Hy((x,R) L) ooyc) is easy to compute; 

H v ((x, R) \Joj V c) = x 2vR + 6 2 -kR 2 + 6> 3 

(it is not worth using inequalities of Lemma f4 . 1 2 1 here) . So, Hy((x, R)L)ujyc) is bounded by 
a positive constant K only depending on parameters 0\, 02, 03 and radius R\. Henceforth, 

P {#^Ux[R ,R +e] = 1) ^V\U = I 

= - i . / e -H v {( X> R)Uu vc ) ze -zX(V) X{dx)Q{dR) 

Zy{uyc) J Ux [ RotRo+e ] 



p -z\(V) 

> rZe- K X{U)Q([Ro,Ro + e]) . 

Ay [LOyc ) 



Recall that Q([Rq, Rq + e]) is positive by (|2,4p . Using the identity (|4.4|l . we finally upper- 
bound the conditional probability P{ooy = 0|wyc) by 

(ze- /< A(C/)Q([ J R ,i? + e]))^ 1 i :, (#wc/x[i?o,/?o+e] = u v \u = $\u V o) . 
This proves Lemma IOI with C = (e~ K X(U)Q([R , R + e]))- 1 . □ 



4.2 Proof of the Connection Lemma 
4.2.1 Outline 

Let us recall that N^:(oo) denotes the number of connected components of Co a having at 
least one ball centered in one of the four cardinal boxes B^, B$, B E or By. In this section, 
we assume N^(lo) > 2. Our strategy consists in exhibiting a subset B of the diamond box 
A in which cob = 0- Moreover, for x £ B, if we are able to control uniformly the energy 
Hb((x,R) U Wflc) on oob c , then the configuration lo should contain a point centered in B 
with large probability as z tends to infinity. This leads to the Connection Lemma. 

For x E B, let us denote by Mhol((x, R),uob c ) the hole number variation when the ball 
B(x,R) is added to the configuration oob c - This quantity is central in our proof. Indeed, 
a first upperbound for the energy Hb((x, R) U oob c ) is given by Lemma 14.121 

H B ((x,R)Uoj B -) = h((x,R),u) B c) < K -0 3 N ho i((x,R),Lo B c) , (4.5) 

where if is a positive constant only depending on parameters 0\, 02, 03 and radii Rq, R\. 
So, to upperbound the energy Hb((x,R) U oob c ) it suffices to upperbound the number of 
created holes (resp. deleted holes) when #3 is negative (resp. positive). This is the reason 
why the proof of the Connection Lemma differs according to the sign of the parameter #3. 
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4.2.2 When #3 is negative 

Let w be a configuration and a be a positive real number. A couple (x, R) £ R 2 x [Rq, R±] 
is said good if all the connected components of the set wa H B(x,R) have an area larger 
than a. These couples are well-named because adding a ball B(x,R) to the configuration 
WAj with a good couple (x,R), does not create too many holes. 

Lemma 4.5. Let (x,R) £l 2 x [Rq,Ri] be a good couple. Then, 

ttR 2 

Af ho i((x,R),L} A ) < ~ ■ 

a 

Proof. The number of created holes when the ball B{x,R) is added to a; a is smaller than 
the number of connected components of the set Oa H B(x,R). Since (x,R) is good, all 
these connected components have an area larger than a. So, there are at most irR 2 /a such 
connected components. □ 

Let us denote by Bad(ixA,aO the following set: 

Bad(wA, a) = {x £ R 2 , 3R <E [Rq, Rq + e], (x, R) is not good } . 

Lemma 4.6. The area of the set Bad(u}/\,a) tends to as a and e tend to 0, uniformly 
on the configuration wa- 

Lemma 14.61 will be proved at the end of this section. Thanks to Lemmas 14.51 and 14. 6| 
we are now able to prove the Connection Lemma. First, we need a family of (small) non- 
overlapping squared boxes whose union covers the convex hull of the boxes B^, B$, B E and 
B v . Precisely, for k > 0, let us consider a subset B of {v + [0, k[ 2 ,v 6 M. 2 } such that for 
any B,B' in B, BOB' is empty, and 

Conv (B N ,B S ,B E ,B V ) C [j B C A . 
The family B is made up of at most c K = k~ 2 A(A) elements. 

The hypothesis A^(w) > 2 ensures the existence of two elements (x\, •) and (x%, •) of u, 
whose centers x\ and X2 are in the union of the four cardinal boxes B^, B$, B E and B v , 
and whose balls B(x±, •) and B(x2, •) belong to two different connected components of Q, 
say respectively C\ and C2. Let [xi,X2] be the segment in M 2 linking x\ with x<i and d be 
the euclidean distance on M 2 . The continuous application 

/ : x G [x\, X2] ' ^ d(x, C\) — d(x, ui \ C\) 

satisfies f{x\) < and f{x2) > 0. So there exists a point x in [xi,^] such that d(x, C\) 
and d(x,ui \ C\) are equal (and positive). Hence, the ball B(x,Rq) does not contain any 
point of to a- Moreover, since x is in the convex hull of the boxes B^, B$, B E and B v , then 
it belongs to one box of the family B, say B. With n < Rq/^/2, the box B is contained in 
B(x,Ro). Consequently, ujb is empty : 

P (N£{oj) > 2 I co A c) < ^ P (w B = I uj A c) . (4.6) 
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For a given box B G B, let us consider the (random) set J7(wa\b) °f points x G B such 
that, for any radius R G [i?o 5 -Ro + £ ]> the couple (x, it!) is good: 

U(u A \ B ) = B\ Bad(uj A \ B , a) . 

Let x G U(ui A \b) and i? G [i2o>-Ro + £ ]- On the one hand, using (|4.5p . #3 < and Lemma 
14.51 we get 

H b ((x,R)Ulo B c) <K-e 3 M , (4.7) 

where M = M(R\,a) denotes the upperbound given by Lemma 14.51 On the other hand, 
Lemma 14.61 implies the area of U(lo A \b) is larger than k 2 /2 for a and e small enough, 
uniformly on the configuration <^a\_b- It follows: 

P {#UBx[R ,Ro+e] = 1 I w # c ) 

= - } - [ e -H B {{ X ,R)^ B o) ze -zX(B) X{dx)Q{dR) 

Z,b{u B c) jBx[Ro,Ro+e] 



ze zk2 



-I e- HB( - {x ' R)ULUBc) X{dx)Q(dR) 

Z B {u B c) Ju(w A \ B )x[Ro,Ro+e] 



£3 iy* 2 



Lb\ub c ) & 

In the previous inequality, replacing e~ ZK " 2 Zb(wb c )~^ with the conditional probability 
P{ub = ®\^b c ), we obtain 

2e K-e 3 M 

P(U! B =$\u B o) < 



zK 2 Q([R ,R + e]) ' 
Finally, the Connection Lemma derives from the above upperbound, (j4.6|) and with 



c 



^Q([R ,Ro + e}) ■ 



In order to prove Lemma 14.61 we have to locate the set B&d(uj A ,a). Lemma 14.71 says 
that the points (x, .) in Ba,d(u A ,a) are at distance around Rq from ui A . We need some 
notations. Let B(x, R) and B(y, •) be two balls satisfying Rq < R < Rq + e and 

< A(B(x,R)nB(y, •)) < ol . 

Then there exists a positive function g(e, a) tending to as a and e tend to 0, such that 

\d(x,B(y,-))-R \<g(e,a) . (4.8) 

The function g is also allowed to depend on radii Rq and R\. The topological boundary 
du) A is composed of a finite number of arcs. Let a be one of them. This arc is generated 
by an element of the configuration uj a , say (y, •). Now, we can define the circular strip 
S g (a) of width 2g(e,a) by 

S 9 (a) = [x G M 2 ; 3y' G a s.t. X = \ + ^ ~ V \ ^) f > ° and 
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Lemma 4.7. The following inclusion holds; 



Bad(u A ,a) C |J S g (a) . (4.9) 
a, arc of 8uia 

Proof. Let us consider a point x in Bad(wA 5 a)- Let R 6 [Rq,Rq + e] such that (x,R) is 
not good. So there exists a connected component of u>a H -B(cc, i?) of area smaller than a. 
The boundary of this connected component through the open ball B(x,R) is composed of 
a finite number of arcs, say a±, . . . , a n . Let a be one of them realizing the minima 

d(x,a) = min d(x,a,i) . 

l<i<n 

Let (y, •) be the element of the configuration wa generating the arc a. Let S(a) be the 
semi-infinite cone centered at y and with arc a (i.e. the union of semi-line [y, y') for 
y' £ a). Then, x necessarily belongs to S(a). Indeed, the opposite situation could lead 
to the existence of another arc a' satisfying d(x,a') < d(x,a). To sum up, x is in the 
semi-infinite cone S(a) and the area of B(x, R) n B(y, ■) is positive and smaller than a. So 
x satisfies ()4.8p and then belongs to S g (a). □ 

Proof of Lemma \4.6\ Let a be an arc of the boundary BOa- Some geometrical considera- 
tions allow to bound the area of the circular strip S g (a) : 

A(S g (a)) < 4 5 (e,a)length(a) , 

where length(a) denotes the length of the arc a. We deduce from this bound and Lemmas 
ED and [2TQ 

a arc of duiA 

< Y length(a) 
a arc of 

< 4g(e,a)£ A ,{Q A ) with A' = A © B(0, R{) 
2A(A' ® B{0,R )) 



< 4 5 (e,a) 



Rf] 



This latter upperbound does not depend on the configuration oo A . So, this ends the proof 
of Lemma 14.61 □ 



4.2.3 When 63 is positive 

In this section, we still assume that N^(uj) is larger than 2. But this time, our aim consists 
in upperbounding the number of deleted holes when the ball B(x,R), x £ B, is added to 
the configuration ujb c - The existence of a suitable set B derives from Lemma [4.81 Its proof 
is rather long and technical, mainly because of the uniformity of p > with respect to the 
configuration oj. 

Lemma 4.8. Assume N£{u)>2. There exist p > (which does not depend on lo) and 
O = 0{uj) € A such that: 

{%) O is in Conv{B N , B s , B E , By) © B(0, §E ); 
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(ii) B(0, pRo) n uj is empty; 

(Hi) B(0, (1 + p)Ro) does not (totally) contain any hole of Co. 

Let us first explain how the Connection Lemma straight derives from Lemma 14.81 As 
in Section I4.2.2[ we need the family B of non-overlapping squared boxes of length side k. 
But here, B is required to cover a little bit more, i.e. 

Conv (J3,,, B s , B e , Bv) © B(0, ^R ) C [j B , (4.10) 

and parameters k and e are chosen small enough so that 

V2K + e<pR (4.11) 

(where p is given by Lemma [4 . 8 j) . Thanks to statement (i) and (|4.10p . the point O belongs 
to a box B & B. Thanks to (ii), (Hi) and (|4.1ip . ujb is empty and Cob c has no hole in 
B := B © B(0, R + e). Hence, 

P (N£(u>) > 2 | U A c) < P (P (WB = ^ I Wflc) l^ flC has no ho i e in B I W A C ) • (4.12) 

Let us pick a box B £ B, a, couple (x, R) £ B x [Rq, Rq + e] and assume that Gjb c has no 
hole in B. Then, no hole is deleted when B(x,R) is added to lob c - So, the hole number 
variation Nhol(( x , R), u >B c ) is nonnegative. Combining with #3 > and (j4.5|) . the energy 
Hb((x, R) U lvb c ) is smaller than K and we finish the proof of the Connection Lemma as 
in Section First, 



ze 1 



P {#u>Bx[R ,R 0+e ] = l|w B <0 ^ ITT re-^K 2 Q([i? ,#o + £]) • 

Zb[UJb c ) 

Thus, replacing e~ ZK<2 Zb(ojb c )~ 1 by the conditional probability P(ujb = 0| w B c )> we get 

P(uo B = 0| w B c) < 



z K 2 Q([i?o,^o + e]) ' 

Finally, the Connection Lemma derives from the above upperbound, (|4,12p and with 



c e K 

C 



K 2 Q([R ,R + e]) ' 
where c K still denotes the number of boxes contained in the family B. 

Now, let us find a point O and a radius p > satisfying the three properties of Lemma 
14.81 A first applicant for the point O can be obtained following the same method as in 
Section T4.2.2i Based on the hypothesis N^(oj) > 2, this method ensures the existence of a 
point O' in the convex hull of the B$, Bs, B-g, B v J s, such that 

d:=d(O',Ci) = d(O',u A \C 1 )>0 (4.13) 

where C\ denotes a connected component of wa counting by N^(uj). Two cases will be 
considered in the following. In the first one- d > \Rq -the connected components of (Da 
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are far away from O' . So do their holes. Then, the choice O = O' is appropriate. In the 
second case- d < \Rq -we exhibit a region close to O' without hole and choose a suitable 
point O inside. About the radius p, it will be proved in the sequel that any positive real 
number such that 

^2 - 1 



{l + pY < 1 + - , 

V7(l + p)~l < 1 



and 



V7 



+ P\ + 



V 



(i - pY 



V7 



+ p 



< 



VE- i 



(4.14) 
(4.15) 

(4.16) 



is suitable. For instance, p = 0.01 satisfies these three conditions. 



Case 1: d > \R Q . 

By construction, O' is in the convex hull of the boxes _B N , Bs, B E , B v and is at distance at 
least Rq + d from any point x in wa- So, it satisfies properties (i) and (ii) of Lemma [4.81 
Now, let us consider a hole T of Assume in a first time that O' does not belong to T. 
By HSH) and Lemma KTS\ 



d(0',T) 2 > (l + ^J Rl>{l + pfR 2 . 



This means that the hole T is outside the ball B(0' , (1 + p)Ro). Now, assume that O' is 
in T. Since O' is equidistant from two connected components of wa then one of them is 
inside the hole T. Hence, T is too large to be totally covered by the ball B{0' , (1 + p)Ro). 
Consequently, O' also satisfies (Hi). 

Case 2: d < ±R Q . 

Let B(x\, R Xl ) be a ball of the connected component C\ on which the distance d(0\C\) 
is reached. Let us consider the point y\ on the segment [O', sci] satisfying B(yi,Ro) is 
included in B(x\, R Xl ) and 

d(O',B(y 1 ,R )) = d(O',B(x 1 ,R Xl )) = d(0', d) = d. 

In the same way, let us consider a point ?/2 such that B(y2, Ro) is included in a) a \ C\ and 

d(O',B(y 2 ,R )) = d(O',Cd A \C 1 ) = d. 

The region without hole, mentioned at the beginning of the current section and which we 
need, is built from points y\ and y 2 . See Figure [2l Let T> be the infinite line passing by y\ 
and yi- Thus, let us consider two infinite lines T>' and T>" parallel to T> and such that 

d(V,V) = d(V",V) = ^j-R 

(say O' and T>' are on the same side of the line T>). We denote by T~L the intersection of 
the convex hull of balls B(yi,Ro) and B(y2,Ro) with the strip delimited by D' and T>". 
On Figure [2J the border of % is drawn in bold. 
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Lemma 4.9. With the previous notations and hypotheses, there is no hole in %. 

Proof of Lemma \4-9\ The closest hole T to the segment [2/1,2/2] is obtained by pressing a 
ball with radius Rq against B{y\,Ro) and B(y2,Ro). If I denotes the distance between 
T and [2/1,2/2] then 21 is the distance between the center of this pressing ball and [2/1,2/2]- 
Pythagoras Theorem gives (2/) 2 + (Rq + h) 2 = (2R ) 2 in which h denotes 

h ■■= 7^(2/1,2/2) - Ro < d . 

1 /7 

In the worst case, h = t^o- Hence, I is always larger than ^j-Ro, which is the distance 
between T> and T>' . To complete the proof, let us add there is no hole in the balls B(yi, Rq) 
and B(y2,Ro) since they are totally covered by (Da- D 

The idea to conclude the proof can be sum up as follows. The region H is sufficiently 
thick to contain strictly more than half of a ball with radius (1 + p)Rq. Hence, the part 
of this ball outside T~L (this is the hatched region on Figure [2]) has a diameter smaller than 
2Rq. Thanks to Lemma I4.15[ it is possible to choose the center O of this ball so that 
B(0, (1 + p)Ro) fl T~L C does not contain any hole. 

Let T>o be the infinite line parallel to T>" , at distance (1 + p)Ro from T>" and on the same 
side as T> of the line T>" . It derives from (|4.15p that the line T>o is trapped between T> and 
V. Let M be the center of the segment [2/1,2/2]- Let us denote by \z\,zi\ the following 
segment: 

[ Zl ,z 2 ] :=B(M,(l-p)R )nV o - 

See Figure [2j We are going to choose the point O on the segment [21, #2]- To do it, some 
geometrical results about the previous construction are needed. They will be proved at 
the end of the section: 

Lemma 4.10. With the previous notations and hypotheses, the following statements hold: 
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(a) for i = 1,2, d(0',«i)<fi2o; 

(6) ©S(0,pi? ) C5(0',flo + d); 

(c) /or 2 = 1,2, d(yi, m) < (y/S - 1 - p)R . 

By convexity and statement (a), any point of the segment [z\,Z2] is at distance from 
O 1 larger than \Rq. Moreover, O' is in the convex hull of the B^, B$, B^, B v 's. Then, any 
point of [z\,Z2] satisfies the property (i) of Lemma [4,81 

By construction of the point O' , the ball B(0',Rq + d) does not contain any point of oj. 
So does the set [z\, Z2] © B(0, pRo) thanks to statement (b). This means that any point of 
the segment [#1,22] satisfies the property (ii) of Lemma [4,81 

Combining statement (c) with i = 1 and Lemma [4, 141 we check there is no hole of a) a \ Ci 
in the ball B(z\, (1 + p)Rq). Let us run the center of a ball with radius (1 + p)Ro along 
the segment [z±,Z2] from z\ to Z2 until that ball meets a hole of a) a \ C\. Two cases can 
be distinguished. 

• This meet does not happen. Then, the ball B(z2, (1 + p)Ro) does not contain any 
hole of cDa \ C\. It does not contain any hole of C\ either thanks to statement (c) 
with i = 2 and Lemma 14.141 In this case, O = Z2 satisfies the property (iii) of 
Lemma 14.81 

• This meet happens: let O be the corresponding center and T be the corresponding 
hole of (Da \ C\. As just before, the ball B(0, (1 + p)Ro) does not still contain any 
hole of Da \ C\. Denote by C the part of this ball outside H: 

c -.= B{o,(i + p)R ) nn c . 

On the one hand, the diameter of C is smaller than 2Rq thanks to (|4.15p . On the 
other hand, C is pressed against the hole T (there is no hole in H); see Figure [2j 
By Lemma 14.151 the holes of the connected component C\ are at distance from T at 
least 2Rq. So they cannot belong to the set C. Therefore, this point O satisfies the 
property (iii) of Lemma 14.81 

Proof of Lemma \4-lU\ The infinite line T> divides B(M,Rq) into two half-balls; let V be 
the one containing the segment [zi,^]- Since 

d(O', yi ) = d(O',y2) = R + d, 

the half-ball V is included in the ball with center O' and radius Rq + d. This inclusion 
admits two consequences. First, the points z\ and Z2 which are in V, are also in the ball 
B(0', R + d). This implies, for i = 1, 2 

d(0', Zi) < ^R , 

i.e. statement (a). Second, the balls B(zi,pRo) which are included in V, are also included 
in B(0' , Ro + d). So does the set [z±, Z2] ©-8(0, pRo) by convexity. Statement (6) is proved. 
It remains to prove statement (c). Let us introduce the orthogonal projection h\ of z\ over 
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the infinite line V (see Figure [2]). Using d(M, Zl ) = (l- p)R , d(hi,Zi) = (1 + p - ^)R 
and d < \Rq, we get 



d(yi,zi) < 



\ 



+ 



(i - pY 



Rq ■ 



Thanks to ()4.16p . statement (c) follows. 



□ 



4.3 Proofs of geometrical lemmas 

Lemma 4.11. Let A be a bounded closed convex set. For any configuration uj, let us 
denote by £a(w) the perimeter of uj viewed through A: 

£a(w) = C(Q n A) - length(dA n u>), 

where length(dA flw) denotes the lentgh of the boundary of A which is inside the set uj. 
Then, 

2A(A(BB(0,R )) 



£a(w) < 



R() 



Proof. The boundary of u) viewed through A corresponds to a finite union of arcs, say 
{&i)i<i<n- For each arc Oj, coming from the ball B(xi,Ri), we consider the circular strip 
S(ai) of width Rq defined by 



S(a,i) 



x 6 



3x' £ d; s.t. 



x' + p(xi — x') with p > 
and d(x, x') < Rq 



Let us notice that the sets (S(ai))i<i< n are disjoint. Indeed, let suppose that there exists 
x G S(cii) n S(aj) for some i ^ j. Without restriction, we can assume that the distance 
between x and a, is smaller than or equal to the distance between x and a,j. Let y be the 
point on m such that this distance is equal to \y — x\. Then, y has to be strictly included 
in the ball B(xj, Rj) which contradicts the fact that y is on the boundary of Q. 
This allows to compare the sum of the areas of (5 , (aj))i<j< n with A(u)): 

n 2 n 

£a(w) = Y^length(ai) < — )^S(ai) 

tin 



i=l 



i=l 



< 



< 



Ho 

2A(A® B{0,Rq)) 

Rq 



□ 



Lemma 4.12. Let A be a bounded subset o/M 2 , to be a configuration on A and (x,R) be 
an element of A x [Rg,Ri]. Let us denote by A((x,R),uj) the area variation when the ball 
B(x,R) is adding to the configuration uj: 



A((x, R),uj)= A((x, R)Uu)- A(uj) . 
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In the same way, we consider the perimeter variation C((x,R),u) and the connected com- 
ponent number variation N cc ((x , R) , u) . The following inequalities hold. 

< A((x,R),lj) < irRf . (4.17) 

2tt(R 1 + R ) 2 



Rn 



< £((x,R),uj) < 2ttRi . (4.18) 



~ 7T { 1 + ^) <Kc((x,R),uj) < 1 . (4.19) 

Proof. Inequalities (|4.17p . upperbounds of (|4.18p and ()4.19p are obvious. The border length 
of cu which is lost when the ball B(x,R) is adding can be interpreted as the perimeter of 
ui viewed through B(x,R) , i.e. as C§^ x R ^(ui). Thanks to Lemma |4,11[ it is smaller than 

2A(B(x,R)®B{0,R )) 2k{R x + R ) 2 
Ro Ro 

This gives the lowerbound of (|4,18p . It remains to lowerbound J\f cc ((x , R) , ui) . For that 
purpose, the number of deleted connected components when B(x,R) is adding to a), is 
smaller than the number of non-overlapping balls which overlap B(x,R). This number is 
at most 

2tt(R 1 + R ) 
2R 

□ 

Lemma 4.13. Let C be a connected component of wa and T be a hole of C. Any point 
i£l 2 such that x ^ C and x ^ T satisfies 

d(x, T) 2 > d(x, Cf + 2d(x, C)R . 

Proof. Let us consider a connected component C, a hole T and a point x satisfying the 
assumptions of the lemma. Let y be a point of the closure of T such that d(x, T) = \x — y\. 
Necessarily, y is on the boundary of two balls B(z, R) and B{z' , R') of C. Since x belongs 
neither to C nor to T, at least one of z — y or z' — y has a nonnegative scalar product with 
x — y. Say z — y. Given \x — z\ and \y — z\ , the distance \x — y \ is minimal when the vectors 
z — y and x — y are orthogonal. Hence, using \x — z\ > d{x,C) + Rq and \y — z\ > Rq, it 
follows from Pythagoras Theorem that 

d(x,T) 2 > (d(x,C) + Ro? ~ R 2 o , 

which concludes the proof. □ 

The following result is a straight consequence of Lemma 14.131 

Lemma 4.14. Let C, C be two connected components o/wa- Let B{x,R) be a ball of C 
and T' be a hole of C which does not contain B(x,R). Then, 

d(x, T') > V3R . 
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Lemma 4.15. Let T and T' be two holes respectively of two connected components C and 
C o/w A . IfT gtT' and T' <£T then 

d(T, T') > 2R . 

Proof. Let T and T 1 be two holes satisfying the assumption of the lemma. We denote by x 
and y two points belonging respectively to the closure of T and T" such that d(T, T') = \x — 
y\. The point x (respectively y) belongs to the boundary of two balls B(z, R) and B(z', R') 
of C (respectively B{w, r) and B{w' , r') of C). An analysis, as in the proof of Lemma [4. 131 
shows that the distance \x — y\ is minimal in the situation where R = R' = r = r' = Rq 
and z, z' , w and w' form a parallelogram with length side 2Rq. Then the points x and y 
are at the middle of two opposite sides and the result follows. □ 
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